Feasibility of precise navigation in high and low latitude regions under scintillation conditions

Scintillation is one of the most challenging problems in Global Navigation Satellite Systems (GNSS) navigation. This phenomenon appears when the radio signal passes through ionospheric irregularities. These irregularities represent rapid changes on the refraction index and, depending on their size, they can produce also diffractive effects affecting the signal amplitude and, eventually producing cycle slips. In this work, we show that the scintillation effects on the GNSS signal are quite different in low and high latitudes. For low latitude receivers, the main effects, from the point of view of precise navigation, are the increase of the carrier phase noise (measured by s¿) and the fade on the signal intensity (measured by S4) that can produce cycle slips in the GNSS signal. With several examples, we show that the detection of these cycle slips is the most challenging problem for precise navigation, in such a way that, if these cycle slips are detected, precise navigation can be achieved in these regions under scintillation conditions. For high-latitude receivers the situation differs. In this region the size of the irregularities is typically larger than the Fresnel length, so the main effects are related with the fast change on the refractive index associated to the fast movement of the irregularities (which can reach velocities up to several km/s). Consequently, the main effect on the GNSS signals is a fast fluctuation of the carrier phase (large s¿), but with a moderate fade in the amplitude (moderate S4). Therefore, as shown through several examples, fluctuations at high-latitude usually do not produce cycle slips, being the effect quite limited on the ionosphere-free combination and, in general, precise navigation can be achieved also during strong scintillation conditions.Postprint (published version

A method for real-time identification and tracking of traveling ionospheric disturbances using ionosonde data: first results

Traveling Ionospheric Disturbances (TIDs) are wave-like propagating irregularities that alter the electron density environment and play an important role spreading radio signals propagating through the ionosphere. A method combining spectral analysis and cross-correlation is applied to time series of ionospheric characteristics (i.e., MUF(3000)F2 or foF2) using data of the networks of ionosondes in Europe and South Africa to estimate the period, amplitude, velocity and direction of propagation of TIDs. The method is verified using synthetic data and is validated through comparison of TID detection results made with independent observational techniques. The method provides near real time capability of detection and tracking of Large-Scale TIDs (LSTIDs), usually associated with auroral activity.Postprint (published version

Report on the design and specifications of the TID algorithms and products

EU H2020 project TechTIDE deriverable, reporting on the design of the adjustments and upgrades required in order to develop the TID identification algorithms. The added value products that will result from the detection methods will be specified and designed. Each detection method developer will design a validation methodology that will lead to the definition of the confidence metrics.TechTIDE project, funded by the European Commission Horizon 2020 research and innovation program [AD-1], will establish a pre-operational system to demonstrate reliability of a set of TID (Travelling Ionospheric Disturbances) detection methodologies to issue warnings of the occurrence of TIDs over the region extended from Europe to South Africa. TechTIDE warning system will estimate the parameters that specify the TID characteristics and the inferred perturbation, with all additional geophysical information to the users to help them assess the risks and to develop mitigation techniques, tailored to their application. This document is TechTIDE D2.1 “Report on the design and specifications of the TID algorithms and products” and it is an output of TechTIDE Task 2.1 (Specifications for the TID algorithms and the resulting products) of the WP2 (TID identification methodologies) which has the final goal to release the basic algorithms for the TID identification and the value-added products for implementation in the TechTIDE warning system. The document presents the design of adjusted and upgraded TID detection codes, the design of the value-added products, and the validation plan. The design of the adjustments and the upgrades of the different methods are based on the initial requirements gathered among potential users affected by TIDs [RD-1]. Some requirements were brought in from ESA Space Situational Awareness Space Weather (SSA SWE) [RD-2] users' requirements. This way, TID algorithms and product outputs will try to adapt to assess ESA SSA SWE Service Network prerequisites.Preprin

Multi-instrumental observations of the 2014 Ursid meteor outburst

The Ursid meteor shower is an annual shower that usually shows little activity. However, its Zenith hourly rate sometimes increases, usually either when its parent comet, 8P/Tuttle, is close to its perihelion or its aphelion. Outbursts when the comet is away from perihelion are not common and outbursts when the comet is close to aphelion are extremely rare. The most likely explanation offered to date is based on the orbital mean motion resonances. The study of the aphelion outburst of 2000 December provided a means of testing that hypothesis. A new aphelion outburst was predicted for 2014 December. The SPanish Meteor Network, in collaboration with the French Fireball Recovery and InterPlanetary Observation Network, set up a campaign to monitor this outburst and eventually retrieve orbital data that expand and confirm previous preliminary results and predictions. Despite unfavourable weather conditions over the south of Europe over the relevant time period, precise trajectories from multistation meteor data recorded over Spain were obtained, as well as orbital and radiant information for four Ursid meteors. The membership of these four meteors to the expected dust trails that were to provoke the outburst is discussed, and we characterize the origin of the outburst in the dust trail produced by the comet in the year AD 1392.Peer reviewe

Report on TID algorithms

This deliverable presents the TID detection algorithms as improved in response to design principles stated in T2.1 and their testing in the lab environment, verification against measurements taken during quiet and disturbed periods of time, benchmarking for their transition to operations, and final validation to the user requirements of accuracy, timeliness, and coverage.TechTIDE project, funded by the European Commission Horizon 2020 research and innovation program [AD-1], will establish a pre-operational system to demonstrate reliability of a set of TID (Travelling Ionospheric Disturbances) detection methodologies to issue warnings of the occurrence of TIDs over the region extending from Europe to South Africa. TechTIDE warning system will estimate the parameters that specify the TID characteristics and the inferred perturbation, with all additional geophysical information to the users to help them assess the risks and to develop mitigation techniques, tailored to their application. This document is TechTIDE D2.2 “Report on the TID algorithms” and it is an output of TechTIDE Task 2.2 (Development of the TID identification algorithms and products) of the WP2 (TID identification methodologies) which has the final goal to release the basic algorithms for the TID identification and to test a first version of the value-added products for implementation in the TechTIDE warning system. The document highlights four aspects of the TID algorithm release process, (1) Developmentbased on the concept, techniques, and algorithms as stated in TechTIDE D2.1, (2) Verification, an internal testing process that ensures algorithm correctness, (3) Benchmarkingneeded to prepare algorithms to transition to operations, and (4) Validation, an external process of ensuring that developed algorithms are compliant with the stated end user expectations.Postprint (published version

Improved characterization and modeling of equatorial plasma depletions

This manuscript presents a method to identify the occurrence of Equatorial Plasma Bubbles (EPBs) with data gathered from receivers of Global Navigation Satellite System (GNSS). This method adapts a previously existing technique to detect Medium Scale Travelling Ionospheric Disturbances (MSTIDs), which focus on the 2nd time derivatives of total electron content estimated from GNSS signals (2DTEC). Results from this tool made possible to develop a comprehensive analysis of the characteristics of EPBs. Analyses of the probability of occurrence, effective time duration, depth of the depletion and total disturbance of the EPBs show their dependence on local time and season of the year at global scale within the latitude belt from 35°N to 35°S for the descending phase of solar cycle 23 and ascending phase of solar cycle 24, 2002–2014. These results made possible to build an EPBs model, bounded with the Solar Flux index, that simulates the probability of the number of EPBs and their characteristics expected for a representative day at given season and local time (LT). The model results provided insight into different important aspects: the maximum occurrence of bubbles take place near the equatorial anomaly crests, asymmetry between hemispheres and preferred longitudes with enhanced EPBs activity. Model output comparisons with independent observations confirmed its soundness.Postprint (published version

pilot ionosonde network for identification of traveling ionospheric disturbances

Travelling Ionospheric Disturbances (TIDs) are the ionospheric signatures of atmospheric gravity waves (AGWs). Their identification and tracking is important because the TIDs affect all services that rely on predictable ionospheric radio wave propagation. Although various techniques have been proposed to measure TID characteristics, their real-time implementation still has several difficulties. In this contribution, we present a new technique, based on the analysis of oblique Digisonde-to-Digisonde (D2D) "skymap" observations, to directly identify TIDs and specify the TID wave parameters based on the measurement of angle-of-arrival, Doppler frequency, and time-of-flight of ionospherically reflected high-frequency (HF) radio pulses. The technique has been implemented for the first time for the Net-TIDE project with data streaming from the network of European Digisonde DPS4D observatories. The performance is demonstrated during a period of moderate auroral activity, assessing its consistency with independent measurements such as data from auroral magnetometers and electron density perturbations from Digisondes and GNSS stations. Given that the different types of measurements used for this assessment were not made at exactly the same time and location, and that there was insufficient coverage in the area between the AGW sources and the measurement lo cations, we can only consider our interpretation as plausible and indicative for the reliability of the extracted TID characteristics. In the framework of the new TechTIDE project (European Commission H2020), a retrospective analysis of the Net-TIDE results in comparison with those extracted from GNSS TEC-based methodologies is currently being attempted, and the results will be the objective of a follow up paper

St. Patrick’s Day 2015 geomagnetic storm analysis based on Real Time Ionosphere Monitoring

A detailed analysis is presented for the days in March, 2015 surrounding St. Patrick’s Day 2015 geomagnetic storm, based on the existing real-time and near real-time ionospheric models (global or regional) within the group, which are mainly based on Global Navigation Satellite Systems (GNSS) and ionosonde data. For this purpose, a variety of ionospheric parameters is considered, including Total Electron Content (TEC), F2 layer critical frequency (foF2), F2 layer peak (hmF2), bottomside halfthickness (B0) and ionospheric disturbance W-index. Also, ionospheric high-frequency perturbations such as Travelling Ionospheric Disturbances (TIDs), scintillations and the impact of solar flares facing the Earth will be presented to derive a clear picture of the ionospheric dynamicsPostprint (published version

Comportamiento característico de la estructura vertical de la ionosfera en condiciones de calma y perturbadas

Aquesta investigació s'ha centrat en profunditzar en el coneixement del comportament de l'estructura vertical de la regió F de la ionosfera, tant en condicions de calma com pertorbades, i en la seva modelització mitjançant funcions analítiques. Les pretensions d'aquesta investigació han estat motivades per les discrepàncies existents entre les prediccions ionosfèriques del gruix i la forma del perfil de densitat de la regió F en condicions de calma i la seva variació característica, i per l'absència d'un model capaç de reproduir la resposta de l'altura del màxim de ionització en condiciones pertorbades. En aquesta investigació s'ha determinat el comportament patró del gruix i la forma del perfil de densitat electrònica de la regió F en condicions de calma (determinats pels paràmetres B0 i B1 del model Internacional de Referència de la Ionosfera, IRI) en un ampli rang de longituds i latituds. Amb això, s'ha desenvolupat un model global per a cada paràmetre mitjançant una formulació analítica simple que simula les variacions temporals d'aquests en condiciones de calma. La simulació d'aquests models millora (en termes de l'error quadràtic mig, RMSE) les prediccions de l'IRI en un 40% per a B0 i en un 20% per a B1. També s'ha caracteritzat la reacció de l'altura del màxim de ionització, hmF2, a latituds mitges i condicions magnèticament pertorbades, i s'ha determinat un comportament sistemàtic d'aquesta pertorbació, ∆hmF2, la morfologia de la qual depèn del camp magnètic interplanetari (IMF), del temps local, de l'estació de l'any i la latitud. Amb això, s'ha desenvolupat un model empíric que simula la pertorbació d'hmF2 resultant durant tempestes geomagnètiques intenses mitjançant funcions analítiques. Aquest model prediu els esdeveniments d'∆hmF2 amb un 86 % d'encert sense generar falses alarmes i amb un RMSE de 40 km respecte els valors experimentals, que és equivalent al rang de variació experimental obtingut en condicions de calma. Finalment, destacar que també han estat objecte d'estudi en aquesta investigació els mecanismes responsables del comportament ionosfèric tant en condiciones de calma com pertorbades i, especialment, el model de tempesta basat en el paper rector de la circulació del vent neutre termosfèric.Esta investigación se ha centrado en profundizar en el conocimiento del comportamiento de la estructura vertical de la región F de la ionosfera, tanto en condiciones de calma como perturbadas, y en su modelado mediante funciones analíticas. Las pretensiones de esta investigación han estado motivadas por las discrepancias existentes entre las predicciones ionosféricas del espesor y la forma del perfil de densidad de la región F en condiciones de calma y su variación característica, y por la ausencia de un modelo capaz de reproducir la respuesta de la altura del máximo de ionización a condiciones perturbadas. En esta investigación se ha determinado el comportamiento patrón del espesor y la forma del perfil de densidad electrónica de la región F en condiciones de calma (determinados por los parámetros B0 y B1 del modelo Internacional de Referencia de la Ionosfera, IRI) en un amplio rango de longitudes y latitudes. Con esto, se ha desarrollado un modelo global para cada parámetro mediante una formulación analítica simple que simula las variaciones temporales de éstos en condiciones de calma. La simulación de estos modelos mejora (en términos del error cuadrático medio, RMSE) las predicciones del IRI en un 40% para B0 y en un 20% para B1. También se ha caracterizado la reacción de la altura del máximo de ionización, hmF2, en latitudes medias y condiciones magnéticamente perturbadas, y se ha determinado un comportamiento sistemático de dicha perturbación, ∆hmF2, cuya morfología depende del campo magnético interplanetario (IMF), del tiempo local, de la estación del año y de la latitud. Con ello, se ha desarrollado un modelo empírico que simula la perturbación en hmF2 resultante durante tormentas geomagnéticas intensas mediante funciones analíticas. Este modelo predice los eventos de ∆hmF2 con un 86% de acierto sin generar falsas alarmas y con un RMSE de 40 km respecto a los valores experimentales, que es equivalente al rango de variación experimental obtenido en condiciones de calma. Finalmente, resaltar que también han sido objeto de estudio en esta investigación los mecanismos responsables del comportamiento ionosférico tanto en condiciones de calma como perturbadas y, especialmente, el modelo de tormenta basado en el papel rector de la circulación del viento neutro termosférico.The main objective of this research is to improve the knowledge on the vertical structure of the ionospheric F region during both, quiet and disturbed conditions, and its modelling by analytical functions. The main motivations of this research were the existing discrepancies between the predictions of the F region electron density profile thickness and shape during quiet conditions and their characteristic variation, and the absence of a model capable to reproduce the electron density peak height response to disturbed conditions. In this research, the pattern behaviour for quiet conditions of the F region electron density profile thickness and shape (determined by the International Reference Ionosphere model (IRI) parameters B0 and B1) was determined in a wide range of longitudes and latitudes. Then, a global model was developed for each parameter using a simple analytical formulation that simulates their temporal variations during quiet conditions. These model simulations improve (in terms of the root mean square error, RMSE) the IRI predictions by 40 % for B0 and by 20 % for B1. The reaction of the electron density peak height, hmF2, at mid latitudes and magnetically disturbed conditions, was also characterized and the systematic behaviour of this disturbance, ∆hmF2, was determined. The morphology of this disturbance depends on the interplanetary magnetic field (IMF), local time, season and latitude. Furthermore, an empirical model was developed to simulate the hmF2 disturbance during intense geomagnetic storms using analytical functions. This model predicts the ∆hmF2 events with a success of 86 % without generating false alarms and with a RMSE of 40 km with respect to the experimental values, which is equivalent to the experimental variation range obtained during quiet conditions. Finally, the mechanisms responsible of the ionospheric behaviour during both, quiet and disturbed conditions, were also studied in this research, specially the storm model based on the leading role of the thermospheric neutral wind circulation

SSC detection using wavelet analysis

For Sudden Storm Commencements, several methods aiming automatic detection of the events have been proposed. Most of them are based in the analysis of digital geomagnetic field data using some morphological aspects of theses events. Here we present a method using a wavelet analysis on the derivative of H component which is related with the steepness of the rising time. Wavelet analysis is very precise tool for time and frequency determination. A comparative of the results from the detection lists is shown